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STABILITY OF A COLOCATED FINITE VOLUME SCHEME 
FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS 

SEBASTIEN ZIMMERMANN 
(Communicated by Jean-Luc Guermond) 

Abstract. We introduce a finite volume scheme for the two-dimensional incom- 
pressible Navier-Stokes equations. We use a triangular mesh. The unknowns 
for the velocity and pressure are both piecewise constant (colocated scheme). 
We use a projection (fractional-step) method to deal with the incompressibility 
constraint. We prove that the differential operators in the Navier-Stokes equa- 
tions and their discrete counterparts share similar properties. In particular, we 
state an inf-sup (Babuska-Brezzi) condition. We infer from it the stability of 
the scheme. 
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1. Introduction 

We consider the flow of an incompressible fluid in a open bounded set C M 2 
during the time interval [0, T}. The velocity field u : £1 x [0, T] — > K 2 and the 
pressure field p : fl x [0, T) — > K satisfy the Navier-Stokes equations 

(1.1) u t - — Au+(u- V)u + Vp = f, 

Re 

(1.2) divu = 0, 
with the boundary and initial condition 

u|ar2 = 0, u| t=0 = u . 

The terms Au and (u • V)u are respectively associated with the physical phenom- 
ena of diffusion and convection. The Reynolds number Re measures the influence 
of convection in the flow. For equations (|1.1|) - (|1.2|) . finite clement and finite dif- 
ference methods are well known and mathematical studies are available (see [10] 
for example). Numerous computations have also been conducted with finite vol- 
ume schemes (e.g. [T3] and [T]). However, in this case, few mathematical results 
are available. Let us cite Eymard and Herbin [7] and Eymard, Latche and 
Herbin [8] . In order to deal with the incompressibility constraint (|1.2p , these works 
use a penalization method. Another way is to use the projection methods which 
have been introduced by Chorin [1] and Temam [TS] . This is the case in Faure 
[S]. In this work, however, the mesh is made of squares, so that the geometry of 
the problem is limited. Therefore, we introduce in what follows a finite volume 
scheme on triangular meshes for equations (|l.ip - (|1.2p . using a projection method. 
An interesting feature of this scheme is that the unknowns for the velocity and 
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pressure are both piecewise constant (colocated scheme). It leads to an economic 
computer storage, and allows an easy generalization of the scheme to the 3D case. 
The layout of the article is the following. We first introduce (section [2]) some no- 
tations and hypotheses on the mesh. We define (section \2.2\i the spaces we use to 
approximate the velocity and pressure. We define also (section 12. 3p the operators 
we use to approximate the differential operators in (|l.ip - (|1.2p . Combining this 
with a projection method, we build the scheme in section [3] In order to provide 
a mathematical analysis for the scheme, we prove in section 0] that the differential 
operators in (|1.I[) - (|1.2[) and their discrete counterparts share similar properties. In 
particular, the discrete operators for the gradient and the divergence are adjoint. 
Also, the discrete gradient operator is a consistent approximation of its continuous 
counterpart. The discrete operator for the convection term is positive, stable and 
consistent. The discrete operator for the divergence satisfies an inf-sup (Babuska- 
Brezzi) condition. From these properties we deduce in section [5] the stability of the 
scheme. 

We conclude with some notations. The spaces (L 2 ,|.|) and (L°°, ||.||oo) are the 
usual Lebesgue spaces and we set L\ = {q £ L 2 ; J n q(x) dx — 0}. Their vectorial 
counterparts are (L 2 , |.|) and (L°°, y^) with L 2 = (L 2 ) 2 and L°° = (L°°) 2 . For 
k £ N*, (H k , || - ||fe) is the usual Sobolev space. Its vectorial counterpart is (H fe , ||.||fc) 
with H fc = (H k ) 2 . For k = 1, the functions of H 1 with a null trace on the 
boundary form the space Hq. Also, we set Vu = (Vmi, Vu2) T if u = (ui, u 2 ) £ H 1 . 
If X C L 2 is a Banach space, we define C(0, T; X) (resp. L 2 (0,T;X)) as the 
set of the applications g : [0, T] — > X such that t — > |g(£)| is continous (resp. 
square integrable). The norms ||.||c(o,X;X) arL d || .||l2(o,t ; x) are defined respectively 

by ||g||c(o,T;X) = sup fe[0 T] |g(i)| and ||g||z,2(o,T;X) = (j^ \g(t)\ 2 ds"j ' . In all 
calculations, C is a generic positive constant, depending only on f2, Uo and f . 

2. Discrete setting 

First, we introduce the spaces and the operators needed to build the scheme. 

2.1. The mesh. Let Tj, be a triangular mesh of fl: £1 = L>KeT h K. For each 
triangle K e 7^, we denote by \K\ its area and £k the set of his edges. If a € £k 7 
riK,a is the unit vector normal to a pointing outward of K . 

The set of edges of the mesh is £h = ^KeT h £K- The length of an edge a £ £h is \a\ 
and its middle point x CT . The set of edges located inside ft (resp. on its boundary) 
is 4"* (resp. £f l xt ): £ h = £ l h nt U £% xt . If a £ £™\ K a and L a are the triangles 
sharing a as an edge. If a £ £f l xt , only the triangle K a inside is defined. 
We denote by x#- the circumcenter of a triangle K. We assume that the measure 
of all interior angles of the triangles of the mesh are below ^, so that xk £ K. If 
cr £ £™* (resp. a £ £ e h xt ) we set d a = d(x KcT ,x L J (resp. d a = d(x CT , x A -„)). We 
define for all edge a £ £h 



The maximum circumradius of the triangles of the mesh is h. We assume ([6] p. 
776) that there exists C > such that 

Vcr££ h , d{x Kal a) > C\a\ and | a\ > Ch. 

It implies that there exists C > such that 

(2.2) Va££ h , r a >C, 
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and for all triangles K G T~h we have (with a G Ek and h^ ^ the matching altitude) 
(2.3) \K\ = i \a\ h K . a > i \a\ d(x K ,x CT ) > C/i 2 . 

Lastly, if K G Th and L G Ti are two triangles sharing the edge a G £ % h nt , we define 

d(x L ,x CT ) 

olk,l = -77 r- 

d(x K ,x L ) 

Let us notice that ax.L G [0, 1] and + = 1. 

2.2. The discrete spaces. We first define 

P = {q G L 2 ; VAT G T h , q\ K is a constant} . P = (Po) 2 - 

For the sake of concision, we set for all G Po (resp. Vh G Po) and all triangle 
K G Th- qi< = qh\x (resp. vx = v/Ja')- Although P ^ H 1 , we define the discrete 
equivalent of a H 1 norm as follows. For all v/j G Po we set 



1/2 



(2.4) ikiu= Yf ^|v^-vkJ 2 + J2 ^I v *j 2 

where t ct is given by (|2.ip . We have [6] a Poincare-like inequality for Po: there 
exists C > such that for all Vh G Po 

(2.5) |v h | <C\\v h \\ h . 
We also have the following inverse inequality. 

Proposition 2.1. There exists a constant C > smc/i that for all v^ G Po 

h\\v h \\ h <C\w h \. 

PROOF. According to (|2.4p 

We deduce from and ([13]) that h 2 r a < C\K a \ and h 2 T a < C\L a \. Thus, 

since \v La - v K J 2 < 2 (|v i<7 1 2 + (v^ | 2 ) , we get 

^llv.H 2 <C Y (\K a \W K J 2 + \L a \\v L J 2 )+C Y I^IIvkJ 2 . 

Hence /i 2 ||v h || 2 < CJ2 KeTh \K\ \v K \ 2 < C | V/l | 2 . 

From the norm \\.\\h we deduce a dual norm. For all v/, G Po we set 

(2-6) ||v ft ||_i,ft= sup . 

For all u/i G Po and G Po we have (u^,,v^) < ||uft||_i,ft ||vft||/i. Now we introduce 
some operators on Po and Po- We define the projection operator IIp : L 2 — * Po 
as follows. For all w G L 2 , llp w G Po is given by 

(2.7) V K G T h , (n Po w)| K = ^-L^ w (x)dx. 
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We easily check that for all w G L 2 and v/, G Po we have (IIp w,v/,) = (w, v^). 
It implies that Ilp is stable for the L norm. We define also the interpolation 
operator Ilp : H 2 — > F - F° r all Q € -ff 2 , IIp (7 6 P is given by 

V-fif eT h , Hp q\K = q{xi<)- 

According to the Sobolev embedding theorem, q G H 2 is a.e. equal to a continuous 
function. Therefore the definition above makes sense. We also set IIp = (IIp ) 2 . 
The operator Hp a (resp. Ilp ) is naturally stable for the L°° (resp. L°°) norm. 
One also checks ([2] and [16]) that there exists C > such that 

(2.8) |v-IIp v|<Cfc||v||i, \q-U Po q\<Ch\\q\\ 2 

for all v G H 1 and q G H 2 . 

We introduce the finite element spaces 

Pf = {d £ L 2 ; VJC £ T h , u| K is affine} , 

Pr = HeP^V^e 4 ni , v h \ K „ (x ct ) = « fc | ia (x CT ) , 

Pi = {v/j G (Pi d ) 2 ; v ft is continuous and v h \on = 0}. 

We have PJ C Hj. We define the projection operator IIpj : Hq — * Pf. For all 
v = 0i, V2) G Hj, Ilpjv = {v\,vl) G Pf is given by 

2 2 
= ^) €P?, J] ( V <> V 4) = J] (Vu;,V4). 

i=l i=l 

The operator IIpc is stable for the H 1 norm and ([2] p. HO) there exists C > 
such that for all v G H 1 

(2.9) |v-n Pf v| < Cft||v||i. 



Let us address now the space P™ c . If G P" c , we have usually Vgj, ^ L 2 . Thus 
we define the operator : P™ c -^ P by setting for all qh G P\ c and all K G fh 

VhQhW = T77T / Vg/jdx. 
The associated norm is given by 

IMIv, = (kl' + iv^i 2 ) 172 . 

We also have a Poincare inequality: there exists C > such that for all qh G P™ c nLg 

(2.10) \q h \<C\V h q h \. 

We define the projection operator Upnc. For all qh G P" c , Dp^^ is given by 

(2.11) v^gi 2 , (n Pr% ,0) = ( % ,0). 

We have the following result. 

Proposition 2.2. If qh G -Po> Tlp nc qh is given by 

v<re (n P .. gi )(x„) = q K .- 
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PROOF. For all edge a £ S^, we define the function i/v £ P™ c by setting 

1 if a = a'. 



Vv(x<7') 



otherwise. 



Let us notice that vanishes outside K a U L a if a £ ££ n * and outside K a if 
a £ £% xt ■ Let a £ £™* . Using a quadrature formula we get 

(U Pr q h ,^) = + l^i) (n Prgh )(x CT ) 

and 

Wh,1pa) = (lK ~^T- +qL 



3 3 

For an edge cr € £,f 4 we have (n^cg^, •0 (T ) = l^i (IIp»cg h )(x £r ) and {q h ^a) = 
qK a ^f^- By plugging these equations into (|2.1ip with cp = tp cr ,we get the result. ■ 
We finally introduce the Raviart-Thomas spaces 

RTq ={vh £ Pf ; Vcr € £a, v/j|a- • n/<, CT is a constant, and v h ■ n\on = 0} , 
RT ={v/j £ RTjJ ; V K £ T h , Va££ K , v h \ K „ ■ n K „ <a = v h \ L „ ■ n^}. 

For all v/j £ RT , K £% L and cr S £a w e set (v^ • n^jj = v/j|at ■ n Kiff . We define 
the operator II RTo : H 1 — > RT . For all v € H 1 , II RTo v £ RT is given by 

(2.12) \/ K £T h . y<r££ Kl (n RTo v ■ n K)ff ) CT = ^- / vda. 

W\ Jo 

One checks [3] that there exists C > such that for all v e H 1 

(2.13) |v-n RTo v| <C/i||v||i. 
The following result will be useful. 

Proposition 2.3. For all v S H 1 such that divv = 0, we have I1 R t v £ Po. 

PROOF. Let v h = II R t v and K £ T^. According to [3] there exists ax £ K 2 and 
&K £ K such that: Vx £ K , v^(x) = a.^ + 6a"X. Thus divv/>|/< = 2 6a- On the 
other hand, according to the divergence formula and (|2.12[) 



= / divv<ix= / v-nd7= / v 7, • n dq = I divv/jdx. 
JK JdK JdK Jk 

Hence &a — and we get: Vx £ K , v^(x) = a^. ■ 

2.3. The discrete operators. The equations (|l.ip ~ (|1.2p use the differential op- 
erators gradient, divergence and laplacian. Using the spaces of section [2~2l we define 
their discrete counterparts. The discrete gradient X7 h : P Q — > P is built using a 
linear interpolation on the edges of the mesh (see [16] for details). This kind of 
construction has also be considered in [5]. We set for all qh £ Po and all K £ Th 

^hAhW = ^2 \cr\ [otK„,L a qK a + Ct La . Ka 9L„) njf,, 

(2-14) + -L J2 W\q K „n K ^. 

We have the following result [15] . 

Proposition 2.4. //^ e L 2 is smc/i i/ioi Vhqii = 0, t/ien c/^ = 0. 
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The discrete divergence operator div^ : P — > Pq is built so that it is adjoint to the 
operator V. (proposition 14. 61 below). We set for all qh £ Pq and all K £ Th 

(2.15) div h v h \ K = jj- ^2 W\ (&l„,k„ v K „ + a Ka)La vl^ ■ n Kt<T . 

The first discrete laplacian A. : Pq — > Pq ensures that the incompressibility con- 
straint (|1.2p is satisfied in a discrete sense ( proposition 13. ip . We set for all qh € Pq 

(2.16) A h q h = div h (V h q h ). 

The second discrete laplacian A. : Po — > Po is the usual operator in finite volume 
schemes [BJ. We set for all v. £ Pq and all K £ Th 

In order to approximate the convection term (u • V)u in (|1.1[) we define a bilinear 
form b/j : Po x Po ^ Po using the well-known upwind scheme ([BJ p. 766). For all 
u/j £ Po, Vh £ Po, and all K £ Th we have 

(2.17) b h (u h ,v h )\ K = — Y \ a \ (W ' n K,a) + vk + (u CT • n K . a y v io j. 

We have set u CT = + C(k„,l„ u l ct an d a+ = max(a,0), a~ = min(a, 0) 

for all a £ M.. Lastly, we define the trilinear form bj, : P x P x P ^ R 2 as 
follows. For all £ P , v. £ P , W/, € Po, we set 

(2.18) b/^UfcjVfcjW/,) = |-K"|wk- bfc(u/,, Vft)| K . 
3. The scheme 

We have defined in section [2] the discretization in space. We now have to define 
a discretization in time, and treat the incompressibility constraint (II. 2p . We use a 
projection method to this end. This kind of method has been introduced by Chorin 
[4] and Temam [15]. The basic idea is the following. The time interval [0, T] is 
split with a time step k: [0,T] = Un=o[*">Wil with JVeN* and t n = nk for all 
n £ {0, . . . , TV}. For all m € {2, . . . , TV}, we compute (see equation (|3.2p below) a 
first velocity field u™ ~ u(i m ) using only equation (jl.lj) . We use a second-order 
BDF scheme for the discretization in time. We then project u™ (see equation (13. 4p 
below) over a subspace of P . We get a a pressure field ~ p(t m ) and a second 
velocity field u™ ~ u(t m ), which fulfills the incompressibilty constraint (|1.2p in a 
discrete sense. The algorithm goes as follows. 

First, for all m £ {0, . . . ,N}, we set f™ = ITp f(t OT ). Since the operator n Po is 
stable for the L 2 -norm we get 

(3-1) - |IIp f(i m )| < \{(trn)\ < l|f|l C(0 , T; L 2 )- 

We start with the initial values 

u° h £ P n RT , uj e P n rt rieF n L%. 

For all n £ {1, . . . , AT}, (u" +1 ,p™ +1 , u^ 1 ) is deduced from ujj) as follows. 

• u^ +1 £ Pq is given by 

o f.n+1 _ a n , „n-l -i _ 

(32) _^ ^ — A, i ur 1 +b, l (2<-u»- 1 ,u»+ 1 )+V^ = fr 1 , 
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• pl +1 £ P[ lc n L\ is the solution of 

(3-3) A h (p n h +1 -p n h ) = ±div h u n h +1 , 

• u^ +1 6 Po is deduced by 

(3-4) < +l =K +l - ^h(p n h +1 -Pi)- 

Existence and unicity of a solution to equation (|3.2[) is classical ([6] for example). 
Let us show that equation (|3.3[) has also a unique solution. Let qh € Po H Lq such 
that A/jg^ = 0. According to proposition 14.61 we have for all G Po 

-(A h q h ,q h ) = -(div h(Vh,qh),qh) = {^hQh, Vh,qh) = \V h q h \ 2 . 

Therefore we have Vhqh = 0. Using proposition 12.41 we get qh = 0. We have 
thus proved the unicity of a solution for equation (|3.3p . It is also the case for the 
associated linear system. It implies that this linear system has indeed a solution. 
Hence it is also the case for equation (|3.3p . Let us now prove that for all to S 
{0, . . . , N}, fulfills (TO)) in a discrete sense. 

Lemma 3.1. Ifvh 6 RTo n Po then div/j v/, = 0. 

PROOF. Let K £%. Since g RT , definition (|2. 15|) reads 

div h v h \ K = — ^ kl (iL„,ic + aK,iJ • n^. 
Since a^- r + ar = 1 we conclude that 



div h v h \ K = X] l°1 v ^ ' n *> = V *" ( TjfT 51 H n *> ) = °- 



E 

Proposition 3.1. For all m € {0, . . . , A} we /iai>e div/,. u™ = 0. 

Proof. For to e {0, 1} we have u£ S P n RT and S P n RT . Applying 
the lemma above we get the result. If to G {2, . . . , N}, we apply the operator divh 
to (|3.3p and compare with (|3.4p . ■ 

4. Properties of the discrete operators 

We prove that the differential operators in ||1.1|) - (|1.2|) and the operators defined 
in section T2.3I share similar properties. 

4.1. Properties of the discrete convective term. We define b : H 1 xH 1 — > L 2 . 

For all u e H 1 and v = (vi, v 2 ) € H 1 we set 

(4.1) b(u,v) = (div(w! u), div(t> 2 u)) . 

We show that the operator b/, is a consistent approximation of b. 

Proposition 4.1. There exists a constant C > such that for all v £ H 2 and all 
u £ H 2 fl Hq satisfying div u = 

l|np b(u,v)-b ft (n RTo u,n Po v)||_i^ < cfc||u|| a |M|i. 

Proof. Let u/j = IIrto u an d v/, = IIp v. According to proposition 12.31 we have 
U/j e P . Let A' G 7/, . According to the divergence formula and (|2.7|) we have 

n Po b(u,v)| K = J v(u-n) dcr. 
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i da 



On the other hand, let us rewrite b h (u h , v h ). Let a S £k H • Setting 

_ f v K si (u/, • njf, f ), > 

one checks that w K (u CT -n A> )++v i(T (u CT -ii K ^)~ = v K ,L a (u CT -n^). By definition 
^ • n^ j(T = a L „,K (uk • n K:CT ) + ajf,L CT (u Lct • n^^) ; since u h € RT we get 
u CT • n giJ = (a L<r ,K + uk,l„ ) (% • n K:(T ) = (u K • njf,,) = (u/, • njf !<T ) CT . Using at 
last ([2T2]) . we deduce from (f27TT|> 

Mu^, v/j)|k = 7^7 V] [ v K ,L„{u-n.K,a)da. 

Thus 

(n Po b(u,v) -b h (u A ,v h ))k = 51 / ( v ~ v a',lJ (u • n) cZct. 
Let t/> ft G Po- We have 

(n Po b(u,v) -b h (u h ,v h ),i/>h) = / ( v - v ^,i J ( u ' n ) 1 

(4-2) = ]T (Vk ct - </>lJ / (v - vj^.O (u • n) da. 

Let a S We want to estimate the integral over a. Since we work in a two- 

dimensional domain, we have the Sobolev injection H 2 C L°°. Thus 

/ (v - v Ka ,L a ) (u • n) da < ||u||j~ / \v-v Ka:L Jda < C\\u\\ 2 / \v-v K<r>L Jda. 

Let us first assume that v £ C 1 . We set 

_ f Xrt ct si (u/, • nK :<T ) CT > 
XK<T,i " 1 X£„ si (u h • n K , CT ) CT < 
If x G cr, we have the following Taylor expansion 

v(x)-V K ^=v(x)-v(x^ iL(7 )= f VV (tX + (1 - t) X Kril J (X - Xjf^iJ (It. 

Jo 

We have |x — x^x„ I < ^- Thus, integrating over a and using the Cauchy-Schwarz 
inequality, we get 

/ |v- \ K „,L a \da< [ [ \Vv(tx+(l-t)x K ^ L „)\ 2 hV~tdtda 

Jo \JaJo 

We then use the change of variable (£,x) — * y = tx+ (1 — t) x.K a ,L a - Let D a be 
the quadrilateral domain given by the endpoints of a, x^ and x^. The domain 
[0, 1] x a becomes DK a ,L a with 

D a n K a si (iifc • tijc,^ > 
D a n L cr si (u/, • < 

For all t G [0, 1] we have /iVi < ht < C d(x KayLt7 , a) t thanks to the hypothesis on 
the mesh. We check easily that d(xK <r ,L <T ,(^)tdtda = dy. Thus we get 

W-VK„,Ljda<Ch(j^ |Vv( y )| 2 dyJ . 
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Since (C 1 ) 2 is dense in H 2 , this estimate still holds for v G H 2 . Plugging this 
estimate into (|4.2I) and using the Cauchy-Schwarz inequality we get 

(n Po b(u, v) - b^LTRToU, n Pn v), ip h ) 

<Cfc||u||HM E \*l„-*k.A E / |Vv(y)| 2 dy 



so that 



(n Po b(u,v) -b^nRToU.np^),^,,) < C/i ||u|| H 2 ||V>/Jlv> ||v||i. Using 
then definition (|2.6|) , we get the result. ■ 

Let us consider now the operator b^. Let u G H 1 and v G L°° nH 1 with divu > 0. 

— i i 2 

Integrating by parts we deduce from (|4.1|) : J n v ■ b(u. v) dx = j Q divudx > 0. 
The discrete operator b^ shares a similar property. 

Proposition 4.2. Let ua G Po such that divh > 0. For all v/j G Po we have 

bh(uh, Vft, Vfe) > 0. 

PROOF. Remember that for all edges a G two triangles K a et L CT share a as 
an edge. We denote by K a the one such that • nx a ,u > 0. Using the algebraic 
identity 2 a (a — b) = a 2 — b 2 + (a — b) 2 we deduce from (|2.18|) 



2bfe(uh,Vh,Vft) = 2 ^ \cr\v K <7- (vx a - vzj (u CT • nfi^ jCT ) 

so that 2bfc(u fc ,Vh,, V/,) > X^^t \a\ (\v K *\ 2 - |v i<7 | 2 ) (u CT • n^^). This sum can 
be written as a sum over the triangles of the mesh. We get 

2b h (u h ,v h ,v h )> ^2 l v ^| 2 E ( Uff ■ n ^.^)- 
Using finally definition (|2 . 1 5[) we get 

\K\\v K \ 2 {div h u h )\ K >0. m 

KeT h 

The following result states that the operator bh is stable for suitable norms. 

Proposition 4.3. There exists a constant C > such that for all v/j G Po,. w/, G 
P , G P satisfying div h u h = 

|bh(u/„ Vh, v ft )| < C |uft| || va.Ha ||vft.||ft.. 

PROOF. For all triangle K 6 % and all edge a G £r- H we have 

(u CT • n^ CT ) + vk + (u CT • n KiCT )~ v Lct = (u CT ■ n^) v# - \(u a ■ n K<a )\ (v Lct - w K ). 

This way, we deduce from (|4. T[) b/j(u/j, va, wj,) = S% + S% with 

= E VA ' ' WK E I "' ( U<T ' nA >) ' 

A'er,, aes K n£' h nt 

S2 = - 2J w ^ ■ E i u<T ■ nA >i ( vl - _ va- )- 
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By writing the sum over the edges as a sum over the triangles we get 

$2 = ~ E |c| |U<T • *>-K,a\ (Vl ct - Vjf) • (w io . - Wjr). 

Using the Cauchy-Schwarz inequality we get 

/ \ 1/2 / x 1 

Since ua G Po we have the inverse inequality [6] /iHu/jHoo < C|u/,|. Using ()2.2|) 
and (|2.4[) we have 

E - v^J 2 < C E ^|vl ct -v^| 2 < Clival 2 

and E CTef *"* |wl ct -w A -J 2 <C\\w h \\ 2 h . Therefore |5 2 | <C|u h | \\v h \\ h \\w h \\ h . On 
the other hand we deduce from definition (|2. 15|) 

Si= E l^l(vA'-w K )(div h u ft )| K = 0. 

Aer h 

By combining the estimates for Si and S*2 we get the result. ■ 
4.2. Properties of the discrete gradient. 

Proposition 4.4. There exists a constant C > smc/i that for all qh € Po ; 
fc|Vh<Zh| <C|g fc |- 

Proof. Using (|2.14p and the Minkowski inequality, we have for all triangle K G Th 

wiv^ki 2 < E ^f(^+C)+ E iSa- 

Let us sum over K &Th- Since |<r| < /i, using (|2.3[) . we get 



IW<§ E E 0*i&+iAriflL)+E E 1*1 

Thus /i 2 |V, l% | 2 < CEx 6 r h 1*1 & < Ckl 2 - 

We now prove that V/j is a consistent approximation of the gradient. 
Proposition 4.5. There exists a constant C > swc/i i/iai /or all q G i? 2 

|n Po (Vg) - Vh(U Po q)\ <Ch \\q\\ 2 . 
Proof. Let K E Th- Using the gradient formula and definition (|2.14|) we get 

\K\ (n Po (Vg) - V h (ti Po q)) = f Vodx- \K\ V h (fl Po q) = E 1 

v •'A J^- A 



A 



where we have set for all edge a G £k H £ 



J A 



(q- (&K,L a ?(xx) + fti^if g(x L<T ))^ n^ CT der 
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and for all edge a € £k H ^ x *: ^ = la \1 ~ ?( x a:)) n io Squaring and using 
(13]) we get 



|ir||(np (Vg)-v ft (np g) 



A' 



1 1 o-efK (j£z£ 



a |2 



Summing over the triangles if G Tj, we get 

2 C 



(4.3) 



np (v<z)-v,(np n g) 



El': 

KeT h <t££ k 



a |2 
K\ ■ 



We must estimate the integral terms I^. Let K E 7^. Let us first assume that 
q G C 2 (fl). Let er S £k n For x g cr we have the following Taylor expansions 

q(x K ) = g(x)+Vg(x)-(x K -x) + / H(g) (tx* + (1 -i)x)(x x - x) • (xjc -x) t <ft , 

Jo 

?(x i<r ) = g(x)+V?(x)-(x z ,-x)+ / H(?)(t3CL ff + (l-t)x)(xL,-x)-(x Lff -x)tdt, 

Jo 

Vg(x) = Vq(x K ) - [ VVg (tx K + (1 - i)x) (x K - x) dt. 
Jo 

Plugging the last expansion into the two others and integrating over a we get 



(4.4) 



(4.5) 



/ (g(xjr) - q) da = \<r\ Vg(x K ) • (x* - x*) ~ A a K + B a K , 
f (g(x L J - g) da = \a\ Vg(x K ) ■ (x L „ - x CT ) - + 

J a 



We have set for T e {K a , L a } 



(4.6) 
(4.7) 



A°, 



VVg (tx T + (1 - *) x ) ( x t ~ x ) d< rfcr i 



<T JO 



H(g) (tx T + (1 - i)x)(x T - x) • (x T - x)tdtda. 
One can bound these terms as in the proof of proposition 14.11 We get 
(4.8) \A° T \ 2 <Ch 2 f |VVg(y)| 2 dy, \B^<Gh A j |H(g)(y)| 2 dy. 

Now, let us multiply (|4.4|) by -a^^n^a, (I4.5[) by — Q-l„,k ^k,c and sum the 
equalities. Since + ock,L b = 1 we have 

-Q-L a ,K \ (<?(Xk) - q)nK,<rd(T - OLK,L, \ (g( x L„) g) n K,a dd 
= J (<1 - (&K a ,L a g(xA>) + OLL„,K a g(x ii0 .))^ 11a> d<7 = 1£. 

On the other hand 

~OLK,L a (Xjf - X CT ) • n^cr - Qi„ ,if (x Lct - X CT ) • IlA'.cr = ~~&K,L a a La ,K (da - dcr) = 0. 

Therefore we get I£ = -a £(r>J c (A J + B&) - a K>La {A" La + B£j n K , a . Using 
estimates (|4.8p we obtain 

\I?<\ 2 <Ch± [ (|H(g)(y)| 2 + |VVg(y)| 2 ) ( iy. 
Jd„ 
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We now consider the case a £ £ A H £/i xt - For x £ cr we have 

g(x K ) = g(x)+V?(x) • (x K - x)+ f K(q)(tx K + (1 - t)x)(x x - x) • (x K - x)tdt. 

Jo 

Multiplying by n A , CT and integrating over cr, we get — 1% = Jg- njf, ff + -B A n^.cr with 
Jj£ = J Vc/(x) • (xk — x) <ix. Since |x A — x| < h if x £ cr, using a trace theorem, 
we have 

l^l<Cft 2 ||V g || L » (CT) <C/i 2 ^ (|V<?(y)| 2 + | VVc/(y)| 2 ) dy) ' . 
By combining this estimate with (|4.8p . we get 

|^| 2 <2|J CT | 2 + 2|S-| 2 < C/i 4 / |H( 9 )(y)| 2 dy 

+ C/i 4 / (|V 9 (y)| 2 + |VVc/(y)| 2 )dy. 

The space C 2 (f2) is dense in H 2 . Therefore the bounds for Ij ( still hold for q £ H 2 . 
Plugging these bounds into (|4.3[) we get the result. ■ 

4.3. Properties of the discrete divergence. The operators divergence and gra- 
dient are adjoint: if q £ H 1 and v £ H 1 with v • n|an = 0, we get (v, Vg) = 
— (c/,divv) by integrating by parts. For and div/j we state 

Proposition 4.6. Ifv h £ P and q h £ P Q we have: (v h , W h q h ) = -(q h ,dxv h v h )- 

PROOF. Using (|2.14|) one checks that (v h ,V h q h ) = Y^KeT Ik (Si + S 2 + S 3 ) with 

Si= ^ aK , L ° VRr ' n *> > S2= X! H aK ' L " VL " ' ni <" CT ' 

ae£ K n£ l h nt o-e^nf'"* 

and <S 3 = Effg^nf'"' M v a ■ n Ai(T . Since oki, + ol^k = 1 we have 
Si = ^ \ a \ ~ a L„,K) vk • n K:J 

Since n^^ „ = — nj(,a, we also have 

5*2= |cr| a^r, V£„ • n L<ri<T = - ^ |<t| a K . La w La ■ n K . a . 

Therefore 

(vfe,VftC/ft)=- X |cr|(ai,if <r Vj t - + aK,L <r Vi ir )-DLir i i ff + ^ |cr| v A - • n A ' : L„- 

Using definition (|2.15[) we get 

(Vfc, VftC/ft) = - X |A'| div/, Vft|if + X X |cr| v A • n A -,L„- 

A'GT h KeT h a££ K 

Since E ff£ £ K M v a • n^,^ = v K ■ J2ae£ K I '' n ^"^ = we obtain finally 

(\h, Vft?ft) = - X c/ A div ft Vh\x = -(qh, div h v h ). m 

K£T h 
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The divergence operator and the spaces Lq, Hj satisfy the following property, 
called inf-sup (or Babuska-Brezzi) condition (see [TU] for example). There exists a 
constant C > such that 

(4.9) inf sup -ML^>C. 

q£Ll\{0} v£Hj\{0) IM|lM 

We will now prove that the operator div^ and the spaces P fl L%, P satisfy an 
analogous property. The proof is based on the following lemma. 

Lemma 4.1. We assume that the mesh is uniform (i.e. the triangles of the mesh 
are equilateral) . Then we have for all qh £ Pq 

^hqh = V/ l (n P ™o^). 

Proof. Since the mesh is uniform we have: Vcr £ olk„,l„ = \- Let K £ Th- 

Using definition (|2.14|) and the gradient formula we get 

n I I 

/ \V h qh - Vfc(npnc gfe )J dx = \ fax* + 9lJ n K ,a 

+ X! n K,T ~ X! j {^P[^qh)^K,ada. 

Since qh £ Pq we deduce from proposition 12.21 



/ U P «cq h da = \a\ (L T pncq, l )(x (T ) 

J a 



a\q Ka ii<J££t* t . 



Plugging this into the equation above, we get Vhqh\x — V ' h{J^p^qh)\K ■ ■ 

Lemma 4.2. We assume that the mesh is uniform. There exists a constant C > 
such that 

W r- D <-, T2 (qhydiVhVh) || 

yq h £P r\L^, sup >Ch\\Up r q h \\ hh . 

v h eP \{o} llVftH/j 

PROOF. If q h = the result is trivial. Let q h £ Pq H Ll\{0}. Let v/, = VhQh € 
Po\{0}. Using proposition 14.61 we have 

~(qh,div h v h ) = (vfo,Vfe<?fc) = \Vhqh\ 2 = \^hQ.h\ \vh\- 

Let xn be the characteristic function of SI. Putting ip — xn m (|2.1ip we get 
H P n.cq h £ Lq. So according to (|2.10p and (|4.ip we have 

|Vh?h| = V h {Tlp r q h ) >C\\Up r q h \\ lh . 
On the other hand, according to proposition 12. II |v/,.| > C h ||vhj|h,. Therefore 

~(q h ,div h v h ) > C h\\Ilpncq h \\ l h \\v h \\ h . ■ 

Proposition 4.7. We assume that the mesh is uniform. There exists a constant 
C > such that for all qh £ Pq H L\ 

sup — i— r — n — ^ > c |n Pr% |. 

v h €P \{0} \\Vh\\h 
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Proof. If q h = the result is clear. Let q h e P n £q\{0}. According to ([479")) 



there exists v G Hq such that 



(4.10) divv = -n P ».« fc and \\v\\i < C \IL Pr q h \. 

We set v/j = IIpjv. We want to estimate — {qh, div/,(IIp v/i)) . Since VhQh € Po 
we deduce from proposition [46] 

-(gfc,div h (IIp Vfc)) = (n Po v fc , Vftgft) = {v hl V h q h ). 

Splitting the last term we get 

(4.11) - divk(IIp v/,)) = (v, V h q h ) - (v - v h , V h q h ). 
One one hand, integrating by parts, we get 

(v, Vhqh) = -(n P nog h ,divv) + X! X! / ( n ^r*i) ( v ' nK, ff )(fe. 
According to (|4. 10|) we have — {Hpneq h) divv) = |IIp«cg ft | 2 . Moreover 

Y Y (Rp?°<lh) {v ■ nK,a) da = Y / (n P »c?h) (v • n^.a) dff 

since v|an = 0. Using [5] p. 269 and (|4.10p we have 



^ / (II P ncg h ) (v • IlK.er) dcr 



< C/lUvldHnpne^ll!^ 

< c/i|n JT .gfc|||n J »e« h ||i ffc . 



So we get 

(4.12) {v,V h q h ) > {\U P ncq h \ - C h\\U Pr q h \\ hh )\Up r q h \. 

On the other hand, using lemma I47T1 and the Cauchy-Schwarz inequality 

|(v-v Aj V h «fc)| = |(v- v h ,V h {H Pr q h ))\ < \v-v h \ |Vh(npj»gfc)|. 
Using ([279]) and (|4~T0)) we get 

|v-v h | = |v-n P cv|<C'/i||v|| 1 <C'/i|n Pr g h |. 

Thus 

|(v - v h , Vh?,,)| < Ch\n Pr q h \ \v h {n Pr q h )\ < C h\W Pr q h \ Hn^^Hi,,,. 

Let us plug this estimate and (14. 12[) into (14. 1 1 [) . We get 

-(gfc,div A (np v fc )) > (|n Pr g h | - Ch\\Up r q h \\ hh ) \n Pr q h \. 

We now introduce the norm We have v/ t = II P jv 6 C H 1 . Thus, using 

[6] p. 776, we get ||II Po v^||/ l < C||v^||i. Since n P j is stable for the H 1 norm, we 
deduce from (|47TU)) 

||v h ||i = ||n P c V ||i < Hvld < C|n P ncg h |. 

Therefore ||np v^||ft < C \Tlpncq h \. Using this inequality in (|4.3p we obtain that 
there exists constants C\ > and C2 > such that 

-(?fc,div A (n Po vfc)) > (Ciinpncghi-Ca/ilin^cghH^) ||n Po v A || fc . 

We deduce from this 

sup JVhidiVhVh) > ^ | H J _ C2 h ||n Pr g h || 1>fe . 
v h eP \{o} \\Vh\\h 
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Let us combine this with lemma |4"?2"1 Since 

Vt > , max (Ct , C\ \U Pr q h \ - C 2 t) > — %- \K Pr q h \ 

we get the result. ■ 

4.4. Properties of the discrete laplacian. We first prove the coercivity of the 
discrete laplacian. 

Proposition 4.8. For all Uh G Po and vu G Po we have 

-(A h u h ,u h ) = \\u h \\\ - (AftUft,v h ) < \\vLh\\h \Wh\\h- 

Proof. Using definition (|2.3p and writing the sum over the triangles as a sum over 
the edges, we have 

-(A/jUft.Vfc) = - ^ V K ■ ( ^ T a (UL a ~ VLk) - ^ T-„Uj(J 

T ° ( V ^- ~ W k) ■ (VL La - U K ) + T a W K -W K . 

tree*"* KeT h je£ K n£™' 

We get the first half of the result by taking v/i = u/j. On the other hand, using the 
Cauchy-Schwarz inequality and the algebraic identity ab + cd < V a 2 + c 2 V b 2 + d 2 , 
we get the second half. ■ 
If v € H 2 , we have |Av| < 1 1 v 1 1 2 . The operator A^ shares a similar property. 

Proposition 4.9. There exists a constant C > such that for all v G H 2 



Ah(n Po v) 



< Cllvllo. 



PROOF. Let v h = n Po v. Let K e T h . According to definition (|2.16p 

(4.13) A h v h \ K = -r^T 51 r * ( v ( x O ~ v ( x *)) ~ i^T XI r ^ v ( x ^)- 

Let us first assume that v = (vi,v 2 ) G (C£°) 2 . Let i G {1, 2}. If a £ £ K C\ and 
x G ct we have the Taylor expansions 

«i( x ij = Wi(x)+Vwi(x)-(x LCT -x)+ / H(ui)(txL <r +(l-i)x)(xi .-x)-(xL <J -x)ttft, 

Jo 

Vi(xji-) = Wi(x)+Vt;i(x)-(xK-x)+ / H(vi)(tx A ' + (l-t)x)(xK-x)-(x7^-x)tdi , 

Jo 

Vui(x) = V« i (x ic ) - / VVv^txK + (l-t)x)(x K -x)dt. 
Jo 

The notation H(«j) refers to the hessian matrix of Uj. Plugging the last expansion 
into the two others and integrating over a, we get 

/ M x lJ - Vi(x)) dx = Vvi(-K K ) ■ (x is - x CT ) - + Bl'l , 

J CT 

/ (u,-(xjf) - «i(x)) c?x = V«,(xjf) • (xk - x CT ) - A^ 1 + 

J a 

The terms A^ z and B^ % are the same as in (I4.6[) and (|4.7p . with Mj instead of g. 
We substract these equations. Since x^ a — = d a n^- j0 . we infer from (|2.1[) 

t ct («<(xij - w,(xk)) = Vvi(x K ) ■ n K ,a + -j- {-A a £l + B^l + A a / - B% 1 ) . 
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Let us consider now the case a € £r- Hf^*. If x G a we have the Taylor expansions 

Vi(x K ) =«i(x)+V«i(x)-(xis— x)+ / H(vi)(txj<-+(l-f)x)(xA:-x)-(xjf-x)idt, 

Jo 

Vui(x) = V«i(xjr) - / Wwi(tx K + (l-t)x)(x iir -x)dt. 
Jo 

Since € Cg° we have i>i(x) = 0. We plug the last expansion into the other and 
integrate over a. Since Xif — x CT = —d a n#-,<r we deduce from (|2.ip 



1 



-r CT Vi(x K ) = Vwi(xA') • n^cr + — (A 7 / - B^) 



Thus we get 

— y 

\K\ L^i 



T a Vi Xjf J 



(x£„)-«i(xjr)) -|— r ^ 

= jlH Vv »( x a') ' E M + E 



I?' 



where we have set for all edge ffe&n 

1 



R' 



ctESk 



and for all edge a <E £ K H £/f f : R % a = (A^ 1 - B% 1 ). Since 2 CTe£K I '! n ^ CT = °' 
setting Rcr = (i?*, i? 2 ), we get 

E T <r ( V ( X L„) — V ( X A")) _ 7^7 E T a v(x K ) = — ^ 

Since the space (Cq 30 ) 2 is dense in H 2 , one checks that this equation still holds for 
v £ H 2 . Using (|4.13p we infer from it 



= E 1*1 



K£T h 



AftV/, 



A' 



A£T h 1 1 crefe 



Using estimates (14.61) and (14. 7|) we obtain 



^ C EE / (|VV^| 2 + |H(^)| 2 ) dx<C||v| 



i—1 <TE£h 

5. Stability of the scheme 

We now use the results of section |4] to prove the stability of the scheme. We 
first show an estimate for the computed velocity (theorem 15. ip . We then state a 
similar result for the increments in time (lemma l5.2p . Using the inf-sup condition 
( proposition ^. 7p . we infer from it some estimates on the pressure (theorem I5.2[) . 

Lemma 5.1. For all m G {0, . . . , N} and n S {0, . . . , N} we have 



,m\2 \~m\2 



, m ~ m 1 2 



o. 



Proof. First, using propositions 13. II and 14.61 we get 
Thus we deduce from (|3.4[) 
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Using the algebraic identity 2 a (a — b) — a 2 — b 2 + (a — b) 2 we get 

2 «\ < - uj») = Kf - Kl 2 + |< - <| 2 = o. 

We introduce the following hypothesis on the initial data. 

(HI) There exists C > such that \u° h \ + |uj,| + k\V h pl\ < C. 

Hypothesis (HI) is fulfilled if we set u° = IIr,To u o and we use a semi-implicit 
Euler scheme to compute u^. We have the following result. 

Theorem 5.1. We assume that the initial values of the scheme fulfill (HI). For 
all to £ {2, . . . , N} we have 



\u™\ 2 + kJ2\K\\l<c 



Proof. Let m £ {2, . . . , N} and 77 £ {1, . . . , m — 1}. Taking the scalar product of 
([372]) with 4fcu£ +1 we get 

3 < +1 - 4 u» + ur 1 A _ 4fc ~ +1 +1 

2k ' ' l y Re V h h ' h ' 

(5.1) +4fcb42<-urSu^\ur 1 )+4MV^,u^ 1 )=4Mf«+ 1 ,u^ 1 ). 

First of all, using lemma 15.11 we get as in [T3] 

4fc( V + i 3ur i -4ug+ur i 

\ h ' 2k 
= K +1 | 2 - |<f + 6 K+ 1 - u» +1 | 2 + |2< +1 - <| 2 - |2< - u"- 1 ! 2 
+ |u»+ 1 -2< + u»- 1 | 2 . 

According to proposition O] we have (A h u™ +1 , u£ +1 ) = #; ||u£ +1 || 2 . Also, 
using lemma 15.11 and (|3.4p , we have 

4fc(V^,u« +1 ) = 4fc(V^,< +1 -< +1 ) 
4 A- 2 

= — (|VK +1 I 2 - IVKI 2 - |VK +1 - vp£| 2 ). 

Multiplying (O by 4fcV ft (^ +1 - p£) and using the Young inequality we get 
4fc 2 



3 



|v(K +i -K)l 2 <3K +1 -u 



h 



According to proposition K2l we have 4fcb/>(2u™ - u™ -1 , u™ +1 , u™ +1 ) > 0. At last 
using the Cauchy-Schwarz inequality, (|2.5p and (|3.1|) we have 

Using the Young inequality we get 

4fe (f^ 1 , u£ +1 ) < 3 fc ||u« +1 ||^ + C k ||f || 2 (0 T;L2) . 
Let us plug these estimates into (15. 1[) . We get 



u 



r 1 ! 2 - ki 2 + |2< +i - <i 2 - 12< - ur 1 ! 2 + k +i - 2u£ + ur 1 ! 2 



4 /c 2 

+3 K +1 - < +1 | 2 + k ||u» +1 || 2 + — (IV,^ 1 ! 2 - |Vrf) < Cfe. 
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Summing from n — 1 to m — 1 we have 

m-l m-l ,,2 

Kf + |2u™ - u™- 1 ! 2 + 3 K +1 < +1 \ 2 + k £ \K + X + -f- ivrfl 2 

<C + 4|ui| 2 + |2u / 1 l -u«| 2 + fc 2 |V^| 2 . 

Using hypothesis (HI) we get the result. ■ 
We now want to estimate the computed pressure. From now on, we make the 
following hypothesis on the data 

feC(0,T;L 2 ), f t e L 2 (0,T;L 2 ) , u eH 2 nHj. divu o = 0. 

For all sequence (g m ) mG N we define the sequence ((5g m ) mS N* by setting Sq m — 
q m — q" 1 ^ 1 for m > 1. We set S = (S) 2 . If the data uo and f fulfill a compatibility 
condition [13] there exists a solution (u,p) to the equations (|1.1|) - (|1.2|) such that 

ueC(0,T;H 2 ), u t € C(0,T;L 2 ), Vp £ C(0, T\ L 2 ). 

We introduce the following hypothesis on the initial values of the scheme: there 
exists a constant C > such that 

(H2) K-uol + illuJt-u^Hoo + l^-p^l^C/i, \ui-u° h \<Ck. 

One checks easily that this hypothesis implies (HI). We have the following result. 

Lemma 5.2. We assume that the initial values of the scheme fulfill (H2). Then- 
there exists a constant C > such that for all m € {1, . . . , N} 

1 

k 



(5.2) !WI<C 



Proof. We prove the result by induction. The result holds for m = 1 thanks to 
hypothesis (H2). Let us consider the case m = 2. We set = u^. Let u^ 1 6 Po 
given by 

2 k ~ ~ 

(5.3) u- 1 = 4 u° - 3 u\ + — A hU i - 2 fe b ft (u° , u\ ) - 2 fc V*p£ - 2 * . 



We substract this equation from equation (|3.4|) written for n = 1. Since 

b h {2 u\ - u° , u 2 ) - b h (u° h ,ui) = b(2 ui - 2 u°, u 2 ) + b ft (u°, Su 2 h ) , 
upon setting = u° — u^ 1 , we get 

h 2k h h - - — A h (Sul) + b h (2 u{ - 2 u° , u 2 ) + b h (ul Sul) = St*. 

Taking the scalar product with 4 k Sul we get 

1 



2 (3 <5£i 2 - 4 Sul + K, Sul) - — (A h (Su 2 h ), Su 



r 'l) 



(5.4) +4 k b h (u° h , Sul, Sul) + 4 k b h (2 uj, - 2 u° , u 2 , <5u 2 ) = 4 fc (<5f 2 , <*u 2 ). 
According to proposition 14.31 we have 

4fc|b h (2u£ - 2u°,u 2 ,«5u 2 )| < Cfc |2u£ -2<| ||u£|U ||«#|U ; 
so that, using hypothesis (H2) 

4fc |b ft (2ui-2u°,u 2 ,<5u 2 )| <Cfc 2 ||u||| h ||5u^|U. 
From the Young inequality and theorem 15. II we deduce 



ik\b h (2u{-ulul,ul-ul)\ < ^\\Sul\\l + Ck s \\ulf h < JL\\Sul\\t + Ck 2 . 



A COLOCATED FINITE VOLUME SCHEME FOR THE NAVIER-STOKES EQUATIONS 19 

On the other hand 

6% = il - fk = n Po f (ta) - n Po f (h) = n Po QT f t ( s ) ds 

Since n Po is stable for the L 2 norm, using the Cauchy-Schwarz inequality, we get 



\K\< J** \tt{s)\ds <Vk ^ \{ t {s)\ 2 d,^ ' <V*||f, 



Thus 



4fc|(<*f 2 ,<Su 2 )| <4k\Sf 2 h \\Su 2 h \ <Ck^ 2 \5ix 2 \. 
So that, using (12. 5| and the Young inequality 

4fc|(«£,*u£)| < Cfc 3 / 2 ||5u^|U < ^ ||<5u 2 || 2 + Cfc 2 . 
The other terms in (|5.4p are dealt with as in the prooof of theroem 15.11 We get 
(5.5) \6ul\ 2 <\6u{f + \2Sul-6u° h \ 2 . 

We know f (|5.2|) for m = 1) that |^u^| 2 < Ck 2 . It remains to estimate the term 
1 2 dul - Su° h \ 2 . According to fO|) 

2 A; 



-<K + — A, U i -2kb h (u° , ujt ) - 2 fc V.pi - 2 * f * ; 



by taking the scalar product with 2 ^u^ — and using the Cauchy-Schwarz in- 
equality we get 



< 2fc( |5U » 



1 



2 k , — |A fc ui| + IV^J + |fi|J |2 K - Kl 



(5.6) + 2*|b(ua,ui,2 5uJt-K)|- 

Let us bound the terms between braces. First, we have 



A h ul = A h (ul - ffp u(t!)) + A h (n Po u(ii)). 
On one hand, according to proposition 14.81 

A^ui-ftp^t!))! 2 = (A h (ui-n Po u(ti)),A h (ui -n Po u(ii)) 

< \\A h (ui - fip u(ti))|| h IK - n Po u(ti)|| fc . 



c 



Applying proposition 12. II we get 



A h (ujt-n Po u(ti)) 



< 



;2 |A h (ui - n Po u(ti))| |u]t - n Po u(t!)|. 

Using the embedding L°° C L 2 we have 

\ul - fi Po u(ti)| = |fi Po « - U(h))\ < ||fip K - "(*l))lloo i 

since Ilp is stable for the L°° norm, we get using hypothesis (H2) 
|ui-fip u(*i)l < Wul-uih^^Ch 2 . 



Therefore 



A h (ui-n Po u(ti)) 



< C. And according to proposition 14.91 

A h (n Po u(ti)) <c||u(ti)|| < c ||u|| C(0,T;H 2 )- 

Hence |A^ujj| < C. Let us now bound the pressure term in (|5.6p . We have 

V/iPi = v fc ($ - fLp p(h)) + (v ft (np p(t!)) - n Po Vp(ti)) + n Po vp(ti) 
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According to proposition 14.41 we have 
Using ([278]) we get 



v h (pl - n Po K*i)) < f H - 5pbP(*i)|. 



Vfc(pi - Spbpfti)) < C ||p(tl)|| 2 < C|b||c(0,T;H=). 

Since Po is stable for the L 2 norm we have |IIp Vp(ii)| < |Vp(ti)| < IHIeto.TiH 1 )- 
Using proposition 14.51 to treat last term we get |VhP^| < C. And according to 

(|3~Tj) and ([572]) for m = 1 we have 4i~ + l f /ll < c - We are left with the term 
|b/i(u^, u\, 2 8u\ — <$u°)| in (|5.6p . We use the following splitting 

Mu^ujj = b h (u^ -IlRToUo.uj,) +b ft (n R T Uo,u^ -fi Po u(*i)) 

+ ^(nRToUo.npgU^i)). 

Let us take the scalar product with 2 8u\ — Su^. We get 

bfc(u£, ujt, 2 «u£ - <5u°) = Si + B 2 + B 3 

with 

Si = b h (u a h - n RTo u , u\, 2 5u\ - Su° h ) , 
B 2 = b fc (n R ToUo,u^ - n Po u(ii),2 Sul - <Su°) , 

and 

B 3 = (b h (n RTo u , fi Po u(ti)),2 du\ - duty . 
Applying propositions 12. II and 14.31 we have 

\Bx\ < £ |u° - n RTo u | \\niWh \2 8u\ - 5u°l- 
According to ([278]) and (|2~T3| we have have 

|u"-n RTo u | = |n Po u°-n RTo u | < |n Po u -u | + |uo-n RTo u | < c/i ||u ||i. 

According to proposition 14.81 and (|2.5[) 



-(AfcuJt.uJtJ^IAfcuiHuil^aiAhuJillluilU; 



since |A^u^| is bounded we get ||u^||/, < C. Hence |£?i| < C|25u^ — Su°\. In a 
similar way, using propositions 12.11 and 14. 3[ we get 

\b 2 \ < £ |n RTo u | \ul - n Fo u(h)\ |2K - 5u°l- 

We have |n RTo u | < |n RTo u - u | + |u | < Ch ||u ||i + |u | < C||u ||i. Using 
moreover ([5]) we get I-B2I < C |2<$uJ^ — <5u°|. Lastly using the following splitting 



b^(n RTo uo,n Po u(ii)) = ^b /l (n RTo u ,n Po u(ii)) - n Po b(u ,u(< 1 )) 

+ n Po b(u ,u(ti)) , 
we have B3 = B31 + B32 with 

B31 = (b h (U RTo u ,U Po u(h)) -n Po b(u , u(ti)), 2Sui-Sul 
B32 = (n Po b(u ,u(ii)),2<K - Suty. 

We have 

|s 3 i| < ||bh(n RTo uo,np u(ti)) -n P Jb(uoMh))\\-iM\ 2 ^h-^h\\h 
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So that, using proposition O \B 31 1 < Ch\\u \\ 2 ||u(*i)|| 2 pSu\ - Su° h \\ h . Using 
proposition 12. II we obtain 

|B 3 i| <C||u || 2 ||u|| 

C(0,T;H 2 ) 

Let us now bound -B32. Using the Cauchy-Schwarz inequality and the stability of 
IIp for the L 2 norm, we have 

\B 32 \ < |n Po b(u ,u(*i))| \2Su{-6u° h \ < |b(u , u(t x )) | \2Sui - Su° h \. 

Integrating by parts, we deduce from (|4.ip 
2 

b(uo,u(*i))| <5^|uo-V«i(*i)| < |uo|||u(ti)||a<C|uo||| u l|c(0,T;H 2 )- 

Thus I-B32I < C |2<$u^ — Su^\. By gathering the estimates for B\, B2, B3 we get 

|b fc (uX,ui,2 5ui-K)| <C. 
Thus we have bounded the right-hand side in (|5.6j) . We infer from it 

\2Sn\-Su%\ <Ck. 

Plugging this estimate into (|5.5p and using (|5.2p for m = 1, we get (|5.2j) for m = 2. 
Let me {3, . . . , TV — 1}. We assume that the induction hypothesis is satisfied up 
to rank n = m — 1. Let us substract equation (|3.2p with the same for n — 1. Since 
the operator b/j is bilinear we get 

^ - Re" ^ ' M 2du ft ~ du ft u /> ) 

+ b„(2 < - ur 1 , <k +1 ) + v fc (<5 P g) = ^r 1 . 

Let us take the scalar product with 4/c£u^ +1 . We get 

( Ma;+1 - 4 /; ; + * ur ' ■«»"¥») - M (*><"*">■««.•") 

+4fcb h (25u£ - aurSufc 1 " 1 .*^ 1 ) + 4fcb fc (2u£ - u£-\5u£+\«5u£ +1 ) 
+4fc(V 7l (<^),<K+ 1 ) = 4fc(<5f£ +1 ,<5< +1 ). 
According to proposition 14.31 we have 

|4fcb h (2*uJ* - 5ur\u£+\5u£ +1 )| < Cfc|25< - (Su^l ||u^ +1 |U 
Using the induction hypothesis we get 

|4fcbfc(25u£-5ur\u^\<K +1 )| <Cfc 2 K+i,j<K +1 IU- 

Using the Young inequality and (|5.ip we infer that 

|4fcb /t (25<-5ur 1 >< +1 »^ +1 )| < ^H^H^ + Cfc 2 . 

The other terms are treated like the case m — 2. We finally obtain ()5.2|) . ■ 

Theorem 5.2. W^e assume that the initial values of the scheme fulfull (H2). There 
exists a constant C > such that for all m € {2, . . . , iV} 

m 

^|n Pr p™| 2 <c. 

n=2 
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Proof. Let m € {2, . . . , TV}. We set n = m — 1. Using the inf-sup condition (|4.7|) 
and proposition 14.61 we get that there exists v/j S Po\{0} such that 

(5.7) c|K|U |n Pr ^ +1 | < -(p«+i )d i Vh v fe ) - (v^ +1 ,v h ). 

Plugging (|3.4[) into (|3.2p we have 



so that 



3< +1 ~4u» 



n-1 



2k 



1 

Re 



^ A.u^ 1 - b„(2 K - ur 1 , + f£ +1 . 



n+1 



V/, 



3< +1 -4u» + <- 
2fc 

,71 — 1 -71 + 1 



1 

Rc 



b,(2<-<-\ur i ,v,)+(fr i ,v,). 

Using the Cauchy-Schwarz inequality, (|2.5p and (|3.ip we have 



3< +1 -4u™ 



u,, 



2k 



-i v/, 



< C 



3< +1 -4u« + <- 



2A- 



and 



Thanks to proposition 14.31 and theorem 15. II we have 



\b h (2ui-ui-\ui+\v h )\ < (2KI + K- 1 !) Hu^iuiKiu^ciiu^iUKiu. 



And according to proposition 14.81 we have ^A/jU^ +1 ,v; 

|3u»+ 1 -4u« + ur 1 | 



(V hP % +1 ,v h )<C + C 
Comparing with (|5.7J) we get 

\Tl Pr pl +1 \<C + C 



< u 



+ K +1 IUKIU 



-7T.+ 1 1 



ft ||v/JU. Thus 



2A- 



|3< +1 



4u" 



.n-l i 



2fc 



u 



Squaring and summing from n = 1 to m — 1 we obtain 



m — 1 



n=l 



4fc 2 



The last term on the right-hand side is bounded, thanks to theorem [5Tj And since 



3< +1 



4u? 



we deduce from lemma 

m — 1 



3« 



n+1 



L ) = 3Su 



n+1 
h 



SuV 



4fc 2 



n=l 



A- 2 
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